An Improved Limit on the Muon Electric Dipole Moment 
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Three independent searches for an electric dipole moment (EDM) of the positive and negative 
muons have been performed, using spin precession data from the muon g — 2 storage ring at 
Brookhaven National Laboratory. Details on the experimental apparatus and the three analyses 
are presented. Since the individual results on the positive and negative muon, as well as the com- 
bined result, = (—0.1 ± 0.9) x 10~^^ e-cm, are all consistent with zero, we set a new muon EDM 
limit, \dn\ < 1.9 x lO"'^® e-cm (95% C.L.). This represents a factor of 5 improvement over the 
previous best limit on the muon EDM. 



PACS numbers: 13.40.Em, 12.15.Lk, 14.60.Ef 

A permanent electric dipole moment (EDM) for a par- 
ticle in a non-degenerate state violates both parity (P) 
and time reversal (T) symmetries. Assuming conserva- 
tion of the combined symmetries CPT, where C refers 
to charge conjugation symmetry, T violation implies CP 
violation. Unlike parity violation, which is maximal in 
weak leptonic decays, CP violation has never been ob- 
served in the leptonic sector. Given the small observed 
levels of CP violation (seen only in the decays of neutral 
kaons and _B-mesons), the standard model (SM) predicts 
that the EDMs of the leptons are so small that their 
detection is well beyond experimental capabilities for 
the foreseeable future. Any non-zero experimental value 
would therefore indicate the existence of new physics. 
This distinguishes leptons from strongly interacting par- 
ticles, which could have a measurable SM EDM if the 
9 parameter in the QGD Lagrangian turned out to be 



sufficiently large. Most models purporting to explain 
the baryon asymmetry of the universe require CP vio- 
lation but the observed level in the i3-meson and 
kaon systems is insufficient, suggesting that there must 
be other, as yet undetected sources of CP violation, which 
could produce non-vanishing EDMs. Consequently, ex- 
perimental searches for EDMs are being widely pursued. 
The muon is of special interest because, again for the 
foreseeable future, it is the only particle, outside the first 
generation, for which a precision EDM measurement is 
feasible [2]. 

The best experimental limits on the EDMs of elemen- 
tary particles or atoms, some of which are listed in Ta- 
ble [H are those for the ^^^Hg atom, the electron, and 
the neutron. The measured values are all consistent with 
zero. These stringent limits, which are likely to improve 
over the next few years, have produced some of the most 



2 



TABLE I: Some experimental limits on EDMs. 



Physical System 


Value, Error (e- cm) 


Reference 


^^^Hg atom 


(0.49 ± 1.29 ± 0.76) x 10"^^ 




electron 


(0.69 ±0.74) X IQ-^'' 


a 


neutron 


(-1.0 ±3.6) X 10-2« 


i 


muon 


(-3.7 ±3.4) X lO"^'' 


i 



significant constraints on extensions to the SMjinany 
of which predict relatively large EDMs 0, S, B [lC|. 
The data on the neutron lead to the current limit llj, 
e < 10-1". 

The current limit on the muon EDM set by the 
last muon g ~ 2 experiment at CERN and also shown in 
Table HI is considerably less stringent. In the SM and 
in some of its extensions, the lepton EDMs scale with 
mass. Scaling the measured electron EDM by the ratio 
of the muon to the electron mass implies a limit for the 
muon of dp = (1.4 ± 1.5) x 10"^^ e-cm. However, the 
scaling which could come from new physics is essentially 
unknown. Indeed, some SM extensionspredict that the 
muon EDM is larger than 10^^'^ e-cm |12l |. 

While the primar y o bjective of the BNL Muon {g — 
2) experiment [H, [13, 111, [H, [O, [Hi was to measure 
the anomalous magnetic dipole moment of the muon, a^, 
this paper describes how, as a secondary measurement, 
a new limit on the electric dipole moment of the muon is 
obtained, a factor of 5 improvement over that achieved by 
the CERN experiment. The improved limit is interesting 
in its own right and helps to resolve one ambiguity in 
interpreting the difference between the theoretical and 
measured values of . 

In general, the electric dipole moment, d, and magnetic 
dipole moment (MDM, fl) are given by 



' i^EDM, where, in the approximation /3 • _B « 0, 



2m 



(1) 



where ij describes the size of the EDM, gq/(2m) is the 
gyromagnetic ratio, q and m are the particle's charge and 
mass, respectively, and S' is a unit vector directed along s, 
the true spin vector. For muons circulating in a storage 
ring, such as those used by both the CERN and BNL 
experiments, the spin vector's precession in the muon 
rest frame (MRF) depends on the interaction of its MDM 
with the magnetic field and of its EDM with the electric 
field. Measurements involving the former interaction can 
provide a determination of the magnetic anomaly, while 
those involving the latter can be used to determine the 
EDM. Note that if CPT symmetry is assumed, the rj 
parameter is the same for positive and negative muons 
while d changes sign. 

In the presence of electric and magnetic fields, in the 
laboratory frame of reference, the rate of spin precession 
relative to the muon momentum direction is given by 
the sum of contributions from the MDM and EDM: Co = 



QuB ± -a 



72 - 1 



fi X 



E 



(2) 



In the approximation f3 • E ^ 0, 



WEDM = ^V7r~ /3 X ^ H 

Zm V c 



2mc 



(3) 



where q is the muon charge, B is the main dipole field and 
E is the field of the focusing electrostatic quadrupoles, 
flp — {g — 2)12 w IQ-"^ is the anomalous magnetic mo- 
ment and F is the Lorentz force. Here, uja, the g — 2 pre- 
cession frequency, includes contributions from both the 
Larmor and Thomas precessions. Both F and therefore 
i^EDM are oriented along the radial direction while the 
magnetic field and therefore tDa are directed vertically. 
A non-zero EDM slightly tips the direction of Qa out of 
the vertical direction and slightly increases the precession 
frequency. 

Under a Lorentz transformation, the laboratory mag- 
netic field leads to a large electric as well as a magnetic 
field in the MRF. In fact, since the effect of the induced 
electric field, from /3 x B, on ujedm is much larger than 
that produced by the focusing quadrupole field, the latter 
is ignored. 

The BNL {g - 2) experiment and the earlier CERN[19|. 
[20I experiment were optimized to measure a^. Both 
experiments stored muons with the "magic" momen- 
tum p « 3.094 GeV/c, corresponding to 7 « 29.3 and 
7T K, 64.4 /is, which gives (7^ — 1)"^ — « 0. With 
this choice of beam momentum, the effect of the focus- 
ing electric field on uja is zero, uja is oriented vertically, 
parallel or anti-parallel to B. 

Both the MDM and EDM measurements rely on the 
detection of positrons or electrons from the three-body 
decays of the muons, /z"*" e'^ + Ve + or /i^ 
e~ +De + i^fi- Positive muons were stored during the 1999 
and 2000 data runs, while negative muons were stored in 
the 2001 data run. (Positrons will be used generically to 
refer to both electrons and positrons in the subsequent 
discussion.) The positron laboratory energies range from 
to 3.1 GeV. All positrons with energies in excess of 1.2 
GeV, which are the ones of interest here, are initially di- 
rected within approximately 40 mrad of the laboratory 
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muon momentum direction. Most positrons have mo- 
menta that are less than those of the muons, and are 
swept by the magnetic field to the inside of the storage 
ring, where they can be intercepted by the detector sys- 
tem. 

A consequence of maximal parity violation in the muon 
decay is a large correlation between the direction of the 
positron momentum and the muon spin in the MRF. For 
^+ (fi^) decay, the positron (electron) momentum is pref- 
erentially directed parallel (anti-parallel) to the spin. 

Transforming to the laboratory frame, the number of 
positrons in a given energy range is modulated accord- 
ing to whether the polarization is along or opposite the 
direction of motion of the muon, following the familiar 
functional form 



N{t) = Noe-^\l + Acos{Lot + 'P)), 



(4) 



where — jtq = A^^ is the dilated muon lifetime. For 
the /i"*" (a*~), N{t) is a maximum when the muon polar- 
ization is parallel (anti-parallel) to the muon momentum. 
The observed value of uj is used to deduce a^, under the 
assumption that the effect of the EDM on the magnitude 
of the spin precession can be neglected. See the discus- 
sion associated with Eq. (fTC|) below. 

For a uniform B field and in the absence of an EDM, 
all muon spins would precess at the same rate, regardless 
of trajectory or momentum, except for very small correc- 
tions for betatron motion and for the effect of the electric 
field on those muons whose momenta are not precisely 
"magic" . The precession vector uj would be anti-parallel 
to the vertical magnetic field or, equivalently, the spin 
vectors would precess in the horizontal plane. To inves- 
tigate the implications of a non-zero EDM, the detector 
acceptance and spin motion are described in more detail. 

The five-parameter equation for N{t), Eq.Q, can be 
rewritten in the more general (differential) form, using 
MRF coordinates 



P(a, 9, (f), t)dadfl = n{a) 



1 + a{a)p ■ s{t) 



dadn, (5) 



where a = {E / Emax)^^^ , the fractional positron decay 
energy in the MRF, p is a unit vector along the decay 
positron's momentum, s is the muon spin and a is the 
decay asymmetry. The tildes refer to MRF coordinates: 
X points in the outward radial direction, y points up- 
ward, and z lies along the muon momentum. The corre- 
sponding polar and azimuthal angles are denoted 9 and 
(j), respectively. The decay asymmetry is given by the 
expression 



a{a) = 



2a -1 



3 - 2a 

Alternatively, in integral form, we have 

Nit) = Noe-^\l + A-s{t)) 



(6) 



(7) 



where s is the precessing polarization vector of the muon 
ensemble, A is the average asymmetry vector, the aver- 
age of ap over detected decays. Of course, the asymmetry 
vector of a data sample depends on the MRF energies and 
angles of accepted positrons. For example, if positrons 
were only accepted if they struck detectors above (be- 
low) the storage ring midplane, then A would acquire a 
positive (negative) vertical component. 

In Eq. ([7]), s{t) describes the spin motion. A includes 
the effects of the weak decay on the electron distribution, 
as well as detector acceptance. Consider the measure- 
ment. We define a local, lab-frame, Cartesian coordinate 
system where x and y are directed radially outward and 
vertically up, respectively and z is longitudinal, along the 
motion of the stored muon beam. For the analysis, 
there is no bias in the spectra of positrons selected, above 
or below the midplane, and therefore Ay « 0. More- 
over, there is little detector bias for positrons with nega- 
tive versus positive radial momentum components. Thus 
Ax w 0. The preferential selection of positrons with high 
laboratory frame energies results in a large value of Az- 
Therefore A « A^z, with z longitudinal and Eq. ([7]) re- 
duces to Eq. (g]). 

For a lab-frame decay electron energy fraction / = 
E/Emax, Az is given by 



Az = 



/-I 



5 + 5/ -4/2' 



with a corresponding (relative) number of muons 



N = 



2(l-/)(5 + 5/-4/2) 
I27r 



(8) 



(9) 



ignoring the effect of detector acceptance. The magnetic 
anomaly is derived from the experimentally deter- 
mined value of bj. The statistical uncertainty on uj in 
Eq. (jl]) is inversely proportional to (A^o^o)^^^ where A^o 
and ^0 refer to the ensemble of accepted events. Graphs 
of differential A (or Az , strictly speaking) , A^ and NA^ 
vs. /, the last a statistical figure of merit (FOM), are 
shown in Fig. [TJ Imposing an energy threshold selects a 
subset of decay positrons with a net average longitudinal 
momentum in the MRF, which leads to the oscillation in 
N{t). 

The asymmetry vector formalism can also be used in 
the context of the EDM analysis. The asymmetry for 
decays selected along any direction u in the transverse 
{x, y) plane is given by 



Au = cos ( 



j(l + 4/)V(l-/)/ 
5(5 -t- 5/ - 4/2) ' 



(10) 



where tanc/) = py/px, with px and Py the transverse mo- 
menta and (j) the azimuthal angle. The corresponding 
asymmetry, A^, and differential statistical figure of merit, 
A^A2 , are shown in Fig. [2l The sensitivity to an EDM is 
greatest over a broad range of lab energies, for 0.2 to 0.7 
of the maximum. 
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FIG. 1: Number TV (solid), g — 2 differential asymmetry 
j4(dashed) and g — 2 statistical figure of merit (FOM, dot- 
ted) NA^ vs. laboratory frame energy fraction. 
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FIG. 2: Differential EDM asymmetry, Au (dashed), and sta- 
tistical figure of merit, NA^ (dotted), vs. laboratory frame 
energy fraction. 



dicular to the total precession vector uj = uja + ujedm, is 
tipped out of the orbit plane by an angle 

(5 = tan"i[7//3/(2a)], (12) 

where a is the anomaly and we see that S is approximately 
proportional to the magnitude of the EDM. The spin 
precession is now given by 

s(t) — (— sj_ cos(5sin (a;i + $) + S|| sin(5)a; 
— (s|| cos 6 + s± sin 5 sin {ojt + $))?/ 
+s±cos{ujt + ^)z, (13) 

that is, the average vertical component of the spin po- 
larization oscillates at angular frequency uj, with an am- 
plitude which is proportional to the EDM. The parallel 
and perpendicular components are defined with respect 
to the new direction of uj. If Syo represents the vertical 
component of the beam's polarization at injection, which 
in the g—2 experiment is very small, the maximum excur- 
sion \sy — Syo\ occurs when the polarization is directed 
either along or opposite to the radial direction. To be 
precise, the vertical EDM oscillation leads the (17 — 2) 
oscillation in N{t) by 90 °. This phase difference, which 
is the same for positive and negative muons, is useful in 
suppressing false EDM signals, as discussed below. It is 
also important to note that the phenomenology of the 
EDM-related spin motion is charge independent - if rj is 
the same for positive and negative muons, the tipping of 
the precession plane and the small change in frequency 
are identical. 

If the asymmetry term in Eq. [5] is given the explicit 
form 

A^y{f, e) cos (1) X + A„j{f,0) sin 4>y + A,{f,0) z, 

(14) 

using lab-frame variables / and 9, the corresponding dot- 
product takes the form 



Concerning the spin motion, with no EDM, the only 
significant torque is that produced by the interaction of 
the magnetic moment with the magnetic field. The spin 
vector precesses in the x-z plane relative to the momen- 
tum vector according to LUa ~ —a^j^{q/m)B . For the 
polarization vector as a function of time is given by 



s{t) = — sj_ sin [ujt + $)i; — s\\y + s±_ cos {ut + <&)£, (11) 

where sy and s_l are, respectively, the magnitudes of the 
components parallel and perpendicular to uj. 

With a non-zero EDM, the precession is no longer 
confined to the horizontal plane. When the torque due 
to the motional electric field acting on an EDM is in- 
cluded, there is an additional spin precession cJedm = 
— {r]/2){q/m)f3 x B, which is directed radially in the stor- 
age ring. The plane of spin precession, which is perpen- 



a(/,0)p(/,0,0)-l(i) = 

A,{f,e){s^cos{ujt + ^))+ 
^xyif, 9) cos0(— cos(5sin(wt + $) + sin(5)-f 
Axy[f, 9) sin0(— S|| cos (5 — s±_ s\n5sin{ujt + $)). 



(15) 



The first term on the right-hand side, which arises from 
the longitudinal component of the asymmetry vector, 
has the cos{ujt + $) time dependence characteristic of 
the anomalous precession. The second (radial asymme- 
try) term, receives a contribution from the EDM, but is 
not readily observable because it cannot be distinguished 
from a small shift in the oscillation phase angle. The 
third (vertical asymmetry) term has the sought-for ver- 
tical oscillation, with an amplitude proportional to the 
EDM. 



As noted above, a non-zero EDM increases the total 
spin precession frequency 

uj — WqVi + tan^ 5. (16) 

While the roughly 3 standard deviation difference be- 
tween the measured and theoretically-predicted values of 
a^i [11] could be the result of new physics such as super- 
symmetry or muon substructure, it could also be caused 
by a shift in the value of produced by a non-zero 
EDM. See Eqs. (O, ^ and dTB]). A sufficiently precise 
measurement of the muon EDM could help resolve this 
ambiguity. 

I. EXPERIMENTAL APPARATUS 

Details on the g — 2 storage ring [2l[, detectors [2^.[23j 
and on the anomalousprecession analysis are presented 
elsewhere [H [H, [H, [ll, [H [ll] . The storage ring mag- 
net is 44 m in circumference with a dipole field of 1.45 T. 
A plan view of the storage ring and detectors is shown 
in Fig. [3l Muons are injected nearly tangent to the cir- 
cumference through the field-free region provided by a 
superconducting infiector. The beam is centered on the 
storage region by means of a small magnetic kick applied 
at 90° downstream from the infiector. Data used in the 
anomalous precession and EDM analyses are collected for 
more than ten muon lifetimes during each injection cycle. 
While the details of beam injection and magnet design 
were rather different, the layout of the CERN III storage 
ring and detectors were very similar to those of E821. 

II. MEASURING THE EDM 

The measurement of the muon EDM requires the de- 
tection of oscillations in py, the vertical momentum of the 
decay positrons, which reflect the oscillations in Sy, the 
vertical polarization of the muon beam. The amplitude 
of the oscillations is proportional to the EDM. A mea- 
surement of the average vertical decay angle vs. time, 
which picks out the third term on the right-hand side of 
Eq. (|15|) . provides the most direct indication of possible 
oscillations in Sy. 

In order to optimize the experimental sensitivity to 
a non-zero EDM, the range of positron energies should 
maximize the FOM for oscillations in Sy (see Fig [2]) 
and maximize sj_. Naturally, the beam polarization is 
optimized for the g — 2 measurement. The EDM FOM 
varies slowly with positron energy - a broad range of en- 
ergies around 1.5 GeV is acceptable. It is important to 
note, however, that the acceptances of the {g — 2) detec- 
tor systems are not optimized for the EDM measurement. 
Those acceptances generally rise with increasing positron 
energy, reach a maximum near 2.8 GeV and then fall to 
the energy endpoint at 3.1 GeV. 

If direct tracking measurements are unavailable (as was 
the case in the CERN III experiment), the observation 
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FIG. 3: A plan view of the g-2 storage ring showing the po- 
sitions of the 24 calorimeters, the traceback chambers and 
other devices. The calorimeters are numbered 1-24 starting 
from the injection point and proceeding clockwise. 



of certain rate oscillations with a sin(a;t -I- $) time de- 
pendence (see Eq. ( [T5|)) can be used instead. In this 
case, one should select data subsets with non-zero average 
positron momentum along y, in order to maximize \Ay\. 
In practice, this means selecting decay positron hits ei- 
ther above or below the storage ring midplane. A further 
refinement is to measure the oscillation phase as a func- 
tion of vertical position on a detector. The sm{u)t + <&) 
time dependence of the vertical oscillation is shifted by 
90 degrees relative to the g — 2 precession {cos{ujt + <!>)) 
and the sin0 term provides a sign flip between the sig- 
nals observed above or below the storage ring midplane. 
In both these approaches, the correlation of detected ver- 
tical position to vertical angle, while strong, is reduced 
somewhat by the range of vertical decay positions. 

The CERN III experiment Q mounted two adjacent 
scintillator paddles on the entrance face (where most 
positrons enter) of one of the calorimeters. One pad- 
dle was mounted above the vertical mid-plane of the 
calorimeter, and the other below. An event was counted 
when a signal from a paddle was registered in coinci- 
dence with a calorimeter signal that exceeded a specified 
energy threshold, typically 1.2 — 1.4 GeV. Any oscilla- 
tion in (jpy) would cause a corresponding oscillation in 
the ratio of the sum and difference of count rates in the 
up {N^) and down {N~ ) paddles. 
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it) = 



N+{t)-N-{t) 
N+{t)+N-{t)' 



(17) 



In the presence of an EDM but in the absence of field 
perturbations which could bend decay electron tracks, 
expressions for TV"*" and N~ can be written in terms of the 
asymmetry vector for each data sample. A sub-sample 
of decay positrons which are detected above the orbit 
plane will have a non-zero vertical asymmetry compo- 
nent, — Ayjj + AzZ. Similarly a sub-sample below 
the orbit plane will have = ~^yy + A^z. Assuming 
the gain and acceptance of the upper and lower detec- 
tors are equal and the storage ring and vertical detector 
midplane are identical, Eqs. ([7]) and with S|| =0 

give 



1 



2 

(1 T AyS±^ sin (5 sin (ujt + $) + ^^s^ cos {ujt + $)). 

(18) 

As expected, the term containing Ay is proportional to 
the EDM, oscillates at angular frequency u and is 90° 
out of phase with the ordinary g — 2 precession given by 
the term containing A^. 

Separating out the contribution from the EDM, the 
expression for becomes 

= i7Voe"^*[lT^EDMsin(a;t-^$)-f A^cos(wt-f-$)], 

(19) 

where Aedm ~ AyS±^ sin 5 is proportional to and = 

Equivalently one can define an angle = 
tan~^ {—Aebm/A^) and an overall asymmetry A = 

K + ^EDM and write 



-At 



1 + Acos (wt + $±*)]. 



(20) 



In the presence of an EDM only (again, without field 
perturbations), can be obtained from either a fit to 
the ratio, Eq. P7)) . 



rit) = 



A-EBM sin (mt + ^) 
1 + Af, cos {cut + $) ' 



(21) 



or from separate fits of Eq. (|20p to the data from the 
top and bottom paddles, with the EDM being inferred 
from the magnitude of angle The latter approach is 
better, since many small spin perturbations, for exam- 
ple from spin resonances, change — N~ in Eq. (PT|) 
without changing ^' in Eq. (|20p . Spin resonances, how- 
ever, are very weak in the g — 2 storage ring. Only high 
longitudinal modes of some non-linear field components 
can oscillate, in the MRF, in resonance with spin preces- 
sion. For such resonances, the original constant part of 
the vertical spin component changes slowly. 



Use of the latter approach led to the CERN result [b^ on 
the combined fi+ and fi- EDMs, = (3.7 ±3.4) x 10^^^ 
e-cm. This is consistent with zero, giving c?^ < 1.05 x 
10~^* e-cm at the 95% confidence level. The overall un- 
certainty is evenly split between statistical and system- 
atic errors. 

The systematic error arises mainly from the uncer- 
tainty in the alignment of the detectors relative to the 
muon beam, which, coupled with the spin precession, 
can produce a false EDM signal. Along the trajectory 
from the muon decay point to the detector, positrons 
emitted with outward radial momentum components will 
travel further than those with inward components, and 
will therefore have more time to spread out in the vertical 
direction. When the muon spin is pointing radially out- 
ward, more positrons are emitted with outward radial 
momentum components than inward, and the width of 
the vertical distribution of positrons at the detector en- 
trance face will be larger than when the spin is pointing 
radially inward. The width of the vertical distribution 
of positrons therefore 'breathes' at the spin precession 
frequency uj/ (27r). 

Ideally, the average distribution of the muons in the 
vertical coordinate is symmetric about ?; = 0, and this 
symmetry is generally refiected in the decay positrons ob- 
served by the detectors. Indeed, if the EDM is zero, if the 
initial vertical component of the spin polarization is zero, 
if the vertical position of the dividing line between scin- 
tillator paddles is aligned with y = 0, and if the paddles 
have the same efficiency of positron detection, then the 
vertical distribution of detected positrons is symmetric 
in y at all times and r(t) = 0. If any of these require- 
ments are not met, then the time average value of r(t), 
(r(i)), will be non-zero in general. A non-zero {r{t)) is 
not, by itself, significant. However, the breathing of the 
vertical width of the beam will now introduce oscillations 
in r(t) at frequency uj/ {2tt), in phase with the oscillations 
that would be produced by a true EDM. In other words, 
it would produce a false EDM signal. 

Corrections for detector inefficiency and misalignment 
can be made, using knowledge of the shape of the ver- 
tical distribution as a function of time, the acceptance 
of the paddles as a function of time, the value of s^o, 
and the measured value of {r{t)). However, in the case 
of the CERN HI experiment, which had only two pad- 
dles mounted on a single calorimeter, detailed measure- 
ments of the vertical distribution could not be made. In 
E821, there is more information on the vertical distribu- 
tion because the detectors have more than two-fold ver- 
tical segmentation. However, the relative acceptances of 
the elements are more complicated because they depend 
not only on the efficiency of the scintillators but also on 
the gain stability and relative efficiency, versus time, of 
the top half of the calorimeter relative to the bottom. 
The time dependence of the vertical distribution on the 
detector face is readily simulated, but not the detector 
response, which must take into account the time-varying 
energy and incident angles of the positrons, as well as 
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FIG. 4: Detail view of a detector station. The calorimeter 
consists of scintillating fiber embedded in a lead-epoxy ma- 
trix, with the fibers being directed radially in the storage 
ring. Five horizontal scintillators, the FSD segments, cover 
the positron entrance face of the calorimeter. Each calorime- 
ter is approximately 23 cm (wide) by 14 cm (high) by 16 cm 
(deep). The PSD (not shown) is placed in front of the FSD. 

geometric variations in gain and efficiency of the scin- 
tillator and calorimeter. Limitations on the corrections 
cause these alignment effects to dominate the systematic 
error. 



III. THE BNL EDM MEASUREMENTS 

In addition to the positron calorimeters, three kinds 
of detectors were used to measure the EDM in the BNL 
experiment: 

• The front scintillation detectors (FSD) are used 
to measure any oscillation in the average vertical 
position of the decay positrons as they enter the 
calorimeters. They consist of an array of five hori- 
zontal scintillator paddles which cover the entrance 
faces of roughly half the calorimeters [2§|. See Figs. 
[3] and m About 10 of the FSDs were instrumented, 
typically those mounted on calorimeters where the 
injection-related background was the least severe. 
The increased segmentation over the CERN ar- 
rangement helps improve knowledge of the vertical 
distribution of the positrons, thereby reducing the 
misalignment error, and permits the more sophis- 
ticated analysis which is described below. 

• The position sensitive detectors [l^l (PSDs) are a 
much more finely segmented version of the FSD, 
with horizontal as well as vertical segmentation. 
They cover the positron entrance faces of calorime- 
ters 13, 14, 15, 18 and 24. 

• The traceback wire chamber system (TWC) con- 
sists of a series of drift chambers mounted in front 



of calorimeter 20, along the path of the incom- 
ing decay positrons. See Fig. O The TWCs 
are used to measure the pitch angle of the decay 
positrons. They also provide information on the 
phase space distribution of the muons in the stor- 
age ring. Tracks measured by the traceback system 
are extrapolated back into the storage volume, to 
the position where the momentum points tangent 
to the storage ring, which is a good approximation 
to the point of decay. 

Another difference between the CERN III and BNL 
EDM measurements concerns the stored beam. The 
functions which describe the time spectra (Eq. ( [7|) , for 
example) are constructed with the assumption that while 
the muons themselves are moving, the spatial distribu- 
tion is static. In fact, the muon beam arrives in a bunch, 
which provides a modulation of the decay signal with pe- 
riod 149.2 ns, the time it takes for the bunch to circle 
through the storage ring. This fast rotation signal, which 
is prominent for roughly the first 70 fis, was filtered out 
in both the CERN III and E821 analyses. However, be- 
cause of our direct muon injection technique, there is 
another very important collective beam motion. For sev- 
eral hundreds of microseconds after injection, the stored 
muon beam exhibits a variety of coherent betatron oscil- 
lations collectively referred to as CBO. Among these is 
an oscillation in the radial beam position (measured at a 
fixed point in the ring) at frequency fcto ~ 465 kHz [isj . 
Descriptions of the vertical distributions of positrons on 
detector faces must account for beam motion. When the 
beam radius is larger than average, the width of the ver- 
tical distribution of positrons at the detector is larger 
because of the longer distance traveled, and conversely, 
when the beam radius is smaller than average, the verti- 
cal width is smaller. 

Combined with a detector-beam misalignment, radial 
beam motion can produce the same sort of vertical os- 
cillation in {py) as the anomalous precession of the spin, 
again proportional to the magnitude of the misalignment. 
However, because the CBO-induced frequency is different 
from that produced by an EDM, this oscillation will not 
be mistaken for an EDM signal, but instead can be used 
to calibrate, and ultimately to correct, the misalignment 
error. 

While the FSD and PSD provide EDM measurements 
patterned on that of the CERN HI experiment, the PSD 
provides another, qualitatively different, approach. The 
phase parameter $ in Eq. varies with y, the detec- 
tor vertical position. As explained earlier, for decays in 
which the polarization points outward, the ensemble of 
decay positrons spread out more in the vertical direction 
than when the polarization points inward. As a con- 
sequence, positrons detected far away from the vertical 
mid-plane {\y\ large) will be more likely to come from 
outward decays than from inward decays. Conversely, 
inward decays will be more likely than outward decays 
when \y\ is small. Therefore $ will depend slightly on y. 
If = 0, the plane of the spin precession is in the hori- 
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FIG. 5: Top and side views (not to scale) of the traceback system. Muons decay to positrons in the storage region. The 
positron then must travel through the thin window in the vacuum chamber scallop, through the traceback chambers and into 
the calorimeter. Horizontally lying straw chambers allow the precise reconstruction the vertical angle of the positron. 



zontal plane, and will be symmetric in y. However, IV. TRACEBACK ANALYSIS OF 1999 AND 

if dp ^ 0, the plane of the spin precession is tipped out 2000 /_t+ DATA SETS 

of the horizontal plane, and <&(?/) will not be symmetric 
in y. 

A. The Traceback Detector 



Finally, although the TWCs were originally designed to 
determine the phase space of the stored muons, for use in 
the anomalous precession analysis, their measurement of 
the average vertical decay angle provides an independent 
measurement of d^, one which is largely immune to the 
detector misalignment problem. In the absence of radial 
magnetic fields, the vertical angle of the track as mea- 
sured in the traceback system, is the same as that at the 
moment of decay. A non-zero EDM would be reflected in 
an oscillation of the vertical component of the positron 
momentum, 90° out of phase with the (g-2) number oscil- 
lation. Further details on the EDM measurements made 
with these three detectors are presented in the following 
sections. 



The traceback detector consists of a set of eight, three- 
layer drift tube planes, designed to measure positron tra- 
jectories along their usual decay path out from the stor- 
age volume into the calorimeters. This detector was in- 
stalled during the 1999 and 2000 running periods at a sin- 
gle position in the ring. By analyzing the positron drift 
time spectrum in a straw, the radius from the anode wire 
at which the particle passed is determined. Tracks are fit 
to these drift circles according to the equations of mo- 
tion of a charged particle in a magnetic field. Due to the 
very inhomogeneous magnetic field in the region of the 
traceback detector (10 T/m), tracks must be integrated 
with a very small step size. The relatively high precision 
of the track reconstruction offers a detailed, time varying 
picture of the decay positron trajectories and the stored 
muon beam. In particular, for the positron's vertical an- 
gle of entrance into the detector, which is examined for an 
EDM oscillation signal, the resolution is approximately 
350 /irad. 
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B. Traceback analysis description and results 



Nent = 4849924 
Chi2/ndf = 414.6/396 
N =1.236e+04±16 
Lifetime = 1.14e+05 ± 7e+03 
W = 2973±12 
Phiq2 = 1.695 ±0.003 
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Nent = 4629682 

Cili2 / ndf = 389.4 / 396 

N =1.158e+04 + 15 

Lifetime = 2.1 e+06 + 2.5e+06 

W = 3724± 11 

phig2 =1.918 ±0.002 
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FIG. 6: Traceback analysis: Plots of the data and fits to 
Eq. Q, modulo the anomalous precession period, for each 
running period, as recorded by the traceback system. (f999 
above, 2000 below) . The g — 2 phase parameter (labeled phig2 
in the fit box) is used in subsequent vertical angle fits. 



where Te is an empirical term which parameterizes both 
the effect of muon decay and that of the recovery of the 
chambers, which are disabled during injection. The phase 
is determined to better than 3 mrad (see Fig. [H]). 

Once the precession phase is established, a plot of aver- 
age vertical angle versus time, modulo the g—2 precession 
period, is fit to the function 

0{t) ^ M + Af, cos{ujt + $) + Aebm sin(wi + *). (23) 

where lo remains fixed as before and $ is set by the pre- 
vious fit. The amplitude of the sine term represents the 
EDM signal. 

The measurement of the vertical angle in the lab frame 
must be converted to a precession plane tilt angle in the 
muon rest frame. (See Eq. (|I2p ). The conversion factor 
is determined through simulation. Several sets of simu- 
lated trajectories were generated and reconstructed, each 
having a different value for precession plane tilt. Shown 
in Fig.[7|is the fit EDM amplitude for various input pre- 
cession plane tilts. A 3 firad amplitude vertical oscilla- 
tion is generated for each milliradian of precession plane 
tilt. 
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Chi2/ndf = 2.013/6 

pO =2.988e-06±9.3e-08 
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As indicated above, a non-zero EDM would generate 
an oscillation in the vertical angle of decay positrons, 
90° out of phase with the g-2 number modulation, N{t). 
While there are several mechanisms which might gener- 
ate oscillations in the measured vertical angle, all are in 
phase with the number ocillation, N{t), and this is also 
measured in the traceback system. By fixing the relative 
phase of the EDM and N{t) oscillations, the false EDM 
signals produced by these effects can be minimized. 

The N{t) spectrum is fit first. To minimize the effect 
of periodic disturbances at other frequencies, which de- 
phase and average away when many time bins are com- 
bined, the data are plotted vs. time, modulo the g-2 pre- 
cession period. The spectrum is fit to Eq. ([^ where 
the precession period T ~ 2tt/u! = 4365.4 ns, is fixed by 
the result of the anomalous precession analysis. 

N{t) = e~*/^= {No + W cos(wt + $)), (22) 



FIG. 7: Traceback EDM fit amplitude vs. simulated input 
precession plane tilts. The fit function is a strict proportion- 
ality: the slope is parameter po with no constant term. A ver- 
tical oscillation with amplitude 3 x 10~^ radians corresponds 
to 1 xlO~^ rad of precession plane tilt. 



The 1999 and 2000 data-sets contain, respectively, ap- 
proximately 4.8 million and 4.6 million well-fit tracks. In 
each data set, approximately 15% of the tracks are back- 
ground: mis-constructed tracks or tracks scattered from 
the upstream vacuum chamber or previous calorimeter 
stations. To reduce the level of mis-constructed tracks to 
the percent level requires cuts which would throw away 
most of the data. The more liberal cuts chosen for the 
final data sample do not induce a false EDM signal. The 
time spectrum of the average vertical angle combined into 
one histogram with a length of the precession period, is 
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FIG. 8: Traceback analysis: Average vertical angle vs. time, 
modulo the anomalous precession period, measured by the 
traceback system for each running period (1999 above, 2000 
below). The data are fit to Eq. (|23p . The amplitude of the 
oscillations, Aeom, is proportional to the EDM. 



then fit with the sum of a sine and cosine, as described 
above. Fits for each data-set are shown in Figs. [S] The 
results are (4.4 ± 6.3) x 10"^ j-g^^j oscillation for the 1999 
data-set and (-4.5 ± 6.2) x 10"^ rad for the 2000 data- 
set. Combining the results, an oscillation amplitude of 
(—0.1 ± 4.4) X 10^^ rad is obtained, where the error is 
statistical. 

Many systematic uncertainties have been studied and 
those relevant to the measurement are listed in Table |IT1 
The "radial field" error refers to the fact that an aver- 
age radial magnetic field around the ring would tilt the 
precession plane in the same way as an EDM. The effect 
of a radial magnetic field on the decay positrons, which 
would change the vertical angle of the tracks, can be ne- 
glected. Another uncertainty comes from the geometry 
of the detector. In making their way to the detectors, 
positrons with an intial outward radial momentum com- 
ponent travel further, on average, than those with an 
initial inward radial momentum component. Combined 



with varying vertical angle acceptance for different decay 
azimuths and an off-center muon distribution, a radial 
beam oscillation can appear in the detector data as a 
vertical oscillation. An error on the number oscillation 
phase or period may result in some mixing between the 
number oscillation signal and the precession plane tilt sig- 
nal. The total systematic error is 0.14 x 10~^^ (e-cm). 
For more details see [IJ]. The negligible systematic er- 
rors indicate that the TWC method should be considered 
in any future attempt to measure the muon EDM. 

Since the systematic errors are negligible, the trace- 
back system's measurement of the EDM for the pos- 
itive muon is determined by the value and (statisti- 
cal) error for the vertical oscillation amplitude alone: 
(— 0.04 ± 1.6) X 10^^^ (e-cm), which corresponds to an 
upper limit of 

\df,+ \ < 3.2 X 10~^^ (e-cm) (95% C.L.). (24) 

V. FSD ANALYSIS OF THE YEAR 2000 
DATA SET 

As described in the introduction, a non-zero EDM 
would result in an oscillation of the mean vertical po- 
sition at the g — 2 frequency but 90° out of phase with 
the number oscillation. Therefore, for each detector, the 
mean position of hits on the FSD, matched to calorimeter 
events with energy 1.4 — 3.2 GeV, is plotted versus time. 
(The center tile is not used in the mean.) An example 
from a single station is shown in Fig. |Vl The plot is fit 
to 



fit) = K+[Sg2Sm{ojt) + Cg2Cos{u}t)] (25) 

_ t 

+ e -CBO X [ScBO sin (cJcBo(i - ^o) + ^cbo) 

+ CcBOCOs(wcBo(<-<o) + *CBo)] 

t 

+ Me . 

The constant term K characterizes the vertical offset be- 
tween the beam and the detector. There are sine and 
cosine terms of frequency lo with the phase fixed such 
that the cosine term is aligned with the number oscil- 
lation. Thus, an EDM signal would appear 90° out of 
phase, in the sine term (S'g2); w is fixed to the frequency 
measured in the g — 2 analysis, to is an empirical term, 
fixed to 100 /isec. 

In addition, there is a sinusoidal term with the fre- 
quency of the coherent betatron oscillation(CBO). Hori- 
zontal focusing causes the beam position near any detec- 
tor to oscillate radially at t^cBO ~ 465 kHz. The oscil- 
lation is prominent at early times, but de-phasing causes 
its amplitude to decay with a lifetime tcbo ~ 100 fis. 
The CBO is clearly evident in a plot of the vertical profile 
width versus time. The frequency ujcbOj lifetime tqbOj 
and phase {^gbo) of the CBO are obtained from fits to 
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TABLE II: Table of systematic errors from the traceback analysis. 



Systematic error 


Vertical 


Precession 


False EDM 




oscillation 


plane tilt 


gener- 




amplitude 


(mrad) 


ated 10 




(/^rad lab) 




(e- cm) 


Radial field 


0.13 


0.04 


0.045 


Acceptance 


0.3 


0.09 


0.1 


coupling 








Horizontal CBO 


0.3 


0.09 


0.1 


Number oscillation 


0.01 


0.003 


0.0034 


phase fit 








Precession period 


0.01 


0.003 


0.0034 


Totals 


0.44 


0.13 


0.14 




o 

1 0.02 

< 

c-0.02 
o 

1-0.04 
'o 

VI 

O-0.06 
-0.08 
-0.1 
-0.12 
-0.14 



Chi2/ndf= 156.6/8 

pO =-0.0868 ±0.0019 




FIG. 9: Sample plot of the mean position on an FSD versus 
time in fill with fit to Eq. (|26|) overlaid. Oscillations at the 
3 - 2 (T « 4365 ns) and CBO (T ^ 2150 ns) frequencies are 
visible. The 3 mm vertical offset was later corrected. (See 
Figs, [island \TT\ below). 



FIG. 10: FSD analysis: 3^2 versus station before the beam 
realignment in the year 2000 data period. The offset from 
zero, indicated by fit parameter po, and inconsistencies be- 
tween detectors indicated by the large of the fit, are due 
to the misalignments between the detectors and the beam. 



the width. Those parameters are then fixed in the fit to 
the average position versus time, Eq. [26] The final term, 
with tm fixed to 60 fis, accounts for slow changes in de- 
tector response and possible pulse pileup effects, when 
two pulses arrive within the detector deadtime. 

Approximately two-thirds of the way through the year 
2000 data run, the stored beam was moved downward 
by 2 mm to improve its alignment with the detectors. 
Fig. [TOj shows the fitted value of the EDM term {Sg2) 
versus station number before the alignment. The figure 
shows significant variations between detectors with an 
average of 0.087 mm, corresponding to an EDM nearly 
as large as the CERN limit. Fig. [TT] is a similar plot 
from data taken after the beam alignment. Although the 
amplitudes of oscillation are much reduced after beam 
alignment, there are still unacceptably large variations 
among detectors. 



As explained above, these inconsistencies arise from 
residual misalignment between the detectors and stored 
beam combined with oscillations in the width of the 
positron vertical profile at the g — 2 frequency. While 
data taken after the beam was realigned show a signif- 
icantly smaller oscillation, if the alignment were perfect 
and the EDM zero, there would be no oscillation at all. 
The amplitude of an oscillation caused by an EDM should 
be consistent from detector to detector. 

The CBO oscillation that was described earlier can be 
used to eliminate the large false EDM signal that results 
from detector misalignment. The CBO causes an oscilla- 
tion in the width of the vertical profile because, for exam- 
ple, decay positrons must travel further to the detectors 
when the muons are at their maximum radial position. 
The oscillation in the mean vertical position at the CBO 
frequency noted in Fig. [Vj results from the misalign- 
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FIG. 11: FSD analysis: Sg2 versus station after the beam 
realignment. The overall offset from zero is smaller when 
the beam is better aligned but the inconsistencies between 
detectors, indicated by a which is only somewhat improved, 
remain. The line represents a best fit constant to the data. 



ment of the detector with the plane of the beam. Since 
the vertical oscillation at the CBO frequency is sensitive 
to the detector misalignment, it can be used to correct 
for the corresponding systematic error in the EDM mea- 
surement. 

A plot of the vertical oscillation amplitude at the .g — 2 
frequency (iS'g2) versus the vertical oscillation amplitude 
at the CBO frequency (^cbo) is shown in Fig. [121 The 18 
points shown are data from the 9 FSDs used in the anal- 
ysis before and after beam realignment. The fit to a line 
produces a good x^, showing, as indicated earlier, that 
the amplitudes of the two oscillations are well correlated. 
The y-intercept, where the CBO oscillation disappears, 
corresponds to an EDM oscillation measurement made 
by a perfectly aligned detector, which thus eliminates a 
large systematic error. 

From simulation, the expected vertical oscillation due 
to an EDM is 8.8 ^m per 10"^^ (e-cm). The Sg2 inter- 
cept, 

Sg2{0) = (1.27 ± 5.88) fim, (26) 

implies an EDM measurement of 

\df,+ 1 < (0.14 ± 0.67) X 10"^^ (e • cm), (27) 

using only the statistical error. A nearly identical re- 
lation between oscillation amplitude and EDM limit is 
obtained in the PSD analysis, presented below. 

Systematic errors dominate the measurement. Any 
source of vertical oscillation at either the g — 2 or CBO 
frequency, in the correct phase, is a potential source of 
systematic error. The largest of the errors is due to the 
tilt of the detectors. A tilt in the detector around the 
beam direction combined with a horizontal oscillation in 



FIG. 12: FSD analysis: the g — 2 sine amplitude vs. CBO 
sine amplitude from fits to the mean vs. time. 



the impact position on the detector face, would result in 
an apparent vertical oscillation. Simulation studies of de- 
cay positron tracks reveal no horizontal oscillation on the 
detector face at the g—2 frequency in the EDM phase, but 
it does indicate a horizontal oscillation at the CBO fre- 
quency of amplitude 0.6 mm. Measurements with a level 
established that the average detector tilt was less than 
8.7 mrad (^ ), implying a systematic error of 6.1 /im on 
Sg2- A more direct measurement of detector tilt, which 
arrived at a much more stringent limit, was made with 
the PSD. (See section |Vll below). 

Another systematic error could arise if the electrostatic 
quadrupoles were themselves tilted with respect to the 
storage ring field. In this small component of 

the nominally radial CBO oscillation frequency would be 
manifest as an oscillation in the vertical mean position of 
the decay positrons, even in a perfectly aligned detector. 
Surveys of the quadrupoles indicate a tilt of less than 2 
mrad; the resulting systematic error is 3.9 /im. 

If the average vertical spin component of the muon is 
not zero then neither is the average vertical angle of decay 
positrons. Since there is a longer average path length for 
positrons emitted when the muon spin is outward than in- 
ward, this would result in an apparent vertical oscillation 
at the detector face, in phase with the expected EDM sig- 
nal. The average vertical angle of positrons approaching 
a detector was measured by the traceback detector. The 
angle was combined with the change in path length for 
positrons from decays when the muon is pointed inward 
versus outward, obtained from simulation, to give an es- 
timate of the vertical oscillation at the g — 2 frequency 
due to muon vertical spin. This effect is also present 
at the CBO frequency since positrons from decays when 
the muons are further out have longer path lengths to 
the detector. This leads to a partial cancellation of the 
systematic error, which is 5.1 /xm, overall. 

A radial magnetic field would defiect decay positrons 
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vertically, causing a similar systematic error. Although 
the average radial field in the ring is known to be less 
than 20 parts per million (ppm) because it affects the 
beam position, the local radial magnetic field felt by the 
decay positrons may be larger than 100 ppm. The size of 
the error can be estimated in an analysis similar to that 
used for the muon vertical spin. Once again, the error 
of 2.5 ^m reflects a partial cancellation from the CBO. 
It should be noted that the average part of the radial 
magnetic field affects the spin similarly to an EDM, thus 
changing the 4" in Eq. (|20p . However, the radial field, 
which is only 20 ppm of the main vertical field, affects the 
spin at a level which is two orders of magnitude smaller 
than our experimental sensitivity to an EDM. The reason 
is that in spin dynamics, the average magnetic field effect 
is almost completely canceled by the average field of the 
focusing electrodes (while in beam dynamics, they cancel 
each other completely). And without an average radial 
magnetic field, the average electric field cannot appear 
in the g — 2 storage ring with its homogeneous vertical 
magnetic field. 

Since each of the top and bottom of the calorimeters 
is read out by a different PMT, there is the potential for 
timing and energy calibration offsets. By applying differ- 
ent calibration constants to the top and bottom of each 
detector, the relative gains of the two halves were cor- 
rected to within 0.5%. The residual error results largely 
from the vertical offset of the beam during the year 2000 
data-taking run. For the 2001 run, where the offset was 
much smaller, the corresponding error is about 0.1%. 
(See section IVip . The timing difference was estimated 
using the cyclotron structure of the data seen early in 
each fill, before the beam de-bunches. The time at which 
the muon bunch passes each detector can be seen as a 
rate oscillation with a 149.2 ns period. By measuring the 
time of the peak of each bunch in the top and bottom of 
the calorimeter the average timing difference was found 
to be less than 0.5 ns. To determine the potential effects 
of these asymmetries, analyses were performed with off- 
sets in the timing and calibration intentionally inserted. 
Based on these analyses, systematic errors due to timing 
and calibration offsets were estimated to be 3.2 fim and 
2.8 /Ltm respectively. 

The sensitivity of the FSD tiles to low energy back- 
scatter from the calorimeters ("albedo") may also cause 
a systematic error. Albedo causes multiple tiles to fire 
at the same time, giving the appearance of two positrons 
hitting the detector when there was only one. Albedo is 
not a problem unless there is a difference in the sensitivity 
to albedo in the top and bottom of the detector. A limit 
on the size of this effect was determined by using several 
methods of dealing with apparent double hits (counting 
both, counting neither, counting one at random). No 
deviations in the results were found larger than 2.0 fim, 
which is taken to be the systematic error. 

There are error bars on both abscissa and ordinate in 
Fig. [T^ Our fits, which use an approximation for the 
latter error, introduce an additional systematic error of 



1 /xm. The effect of FSD tile inefficiency and dead time 
was also investigated but the resulting systematic errors 
were less than 1 /xm and so were ignored. 

The systematic effects, given in Table IIIIl are uncor- 
related; the uncertainties are added in quadrature. The 
total error is 11.9 /im, which gives an EDM measurement 
of 

= (-0.1±1.4) X 10"^^ (e-cm). (28) 

and a limit 

\d^+ \ < 2.9 X 10"^^ (e-cm) (95% C.L.). (29) 



VI. PSD ANALYSIS OF 2001 (^i ) DATA 

The Position Sensitive Detector (PSD) [l^l is a finely 
segmented two-dimensional scintillator hodoscope cover- 
ing the positron entrance faces of five detector stations. 
Each hodoscope consists of one plane of 20 scintillator 
tiles directed horizontally and one plane of 32 tiles di- 
rected vertically. Each tile is 8 mm wide and 7 mm thick. 
The PSDs were mounted on the face of the calorimeters, 
providing vertical and horizontal position data. As indi- 
cated in the introduction, two EDM searches were made 
with the PSD, one modeled on that made with the FSD 
and another in which the symmetry of <i> with y was ex- 
amined. 

While the systematic concerns of the first search are 
closely related to those of the FSD, the latter search re- 
quires that the correct t = he established for every 
tile of the detector. Using coincident timing information 
from the calorimeter, the time offsets for every PSD stick 
can be determined, on a run-by-run basis, and the hits in 
the tiles aligned with the calorimeter within a 25 ns time 
window. Moreover, since the phase is a strong function of 
the decay positron energy, an elaborate program of gain 
balancing was also required. First, the overall gain of 
each calorimeter was established, on a run-by-run basis. 
To this end, the very linear region of the energy distribu- 
tion, from about 60 to 90% of the maximum, was fit to 
a straight line. The x-intercept was taken as the energy 
endpoint, 3.1 GeV. To minimize pulse pileup, only data 
taken more than 250 /xsec after injection were included 
in the calibration energy spectra. After energy scale ad- 
justments and corrections for timing offsets, calorimeter 
times and energies could be matched with a PSD vertical 
coordinate. A similar calibration procedure then estab- 
lished position dependent energy endpoints. In this case, 
where because of shower leakage energy endpoints are not 
expected to sit at 3.1 GeV, the spectra were matched to 
the predictions of a complete tracking and detector sim- 
ulation. 

Time-dependent detector response can be misinter- 
preted as an EDM signal. Variations in gain with time 
were studied separately and eliminated by applying a 
time-dependent gain factor, /(t, Y) for each Y tile. The 
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TABLE III: FSD analysis: Error Table for Sg2, the EDM-sensitive parameter. The conversion to error on the EDM is 0.114x 10 

e-cm/nm. 



Effect 


Error (/xm) 


Detector Tilt 


6.1 


Vertical Spin 


5.1 


Quadrupolc Tilt 


3.9 


Timing Offset 


3.2 


Energy Calibration 


2.8 


Radial Magnetic Field 


2.5 


Albedo and Doubles 


2.0 


Fitting Method 


1.0 


Total Systematic 


10.4 


Statistical 


5.9 


Total Uncertainty 


11.9 
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correction function took the form 

fit, r) - (1 + Aie-*/^^)(l + ^3e-(*-^^^'/^'), (30) 

where parameters Ai — depend on the PSD vertical 
coordinate Y. Pulse pileup was corrected, on an average 
basis, in the time spectra. [18] Finally, time dependence 
arising from the cyclotron and vertical betatron motions 
of the beam were minimized with a simple digital filter. 

In order to study systematic errors related to the CBO, 
in 2001 the storage ring was run at several different field 
indices n. Also, approximately two thirds of the data 
were collected when the average radial magnetic field was 
about 10 ppm of the main vertical B field, which pro- 
duced an average vertical beam displacement of 0.6 mm. 
After the radial field was zeroed, the alignment of the 
beam with respect to the detectors was much improved. 
Data taken before and after the magnetic radial field cor- 
rection were analyzed separately. Data were also divided 
into subsets determined by different run conditions. 

A few cuts were made to the data sample of PSD hits 
and calorimeter pulses. For single PSD cluster events 
(no sign of pileup) the energy was restricted to fall be- 
tween 1.4 and 3.4 GeV. Fits to time spectra were started 
no earlier than 32 /xsec after injection. Table IIVI lists 
the total number of events (after cuts and before pileup 
subtraction) for the PSD analyses. 

A. Details of the two PSD Analyses 

1. Phase Method 




-60 -40 -20 20 40 60 

Y [mm] 



FIG. 13: PSD analysis: PSD 15 before change. Phase of 
g — 2 oscillation vs. vertical position, with fit to Eq. (|31|l . 
pi, the EDM-sensitive parameter, describes any asymmetry 
between the two sides of the graph.. 



serves as a cross-check on the phase vs. vertical posi- 
tion analysis. For each PSD station, the average vertical 
position and rms width were fit vs. time. The rms ver- 
tical width was fit first in order to determine the CBO 
period and lifetime parameters (Tcbo and Tcbo), which 
were then fixed in the fit of the mean vertical position. 
The rms vertical width distribution vs. time, was fit to 
the 9-parameter function: 



The ratio method [18l| was used to determine the 
g — 2 phase vs. vertical coordinate for each of the 20 
horizontally-directed tiles on each detector. The three 
parameter (A, to and $) fits of the ratio method are par- 
ticularly insensitive to any residual slow variations in de- 
tector performance. The phase vs. vertical coordinate 
plot was then fit to the 4-parameter function 



'^{y) = Po + Pi{y - P2) + \Pz{v-P2)\ 



(31) 



where pi, which is proportional to the EDM, describes 
the up-down asymmetry (see [2^ for details including 
the determination of the proportionality constant used 
below). Parameter p3 accounts for phase changes not 
related to the EDM signal, p2 corresponds to the vertical 
detector offset and po is an overall detector phase offset. 
The result of the fit for one detector, for one data set, 
is plotted in Fig. [T3] and the pi-variable for all PSD fits, 
both before and after the radial magnetic field correction 
are plotted versus PSD station in Fig. [HI 



2. Vertical Position vs. Time Method 

The vertical position vs. time method is modeled 
closely on that used in the analysis of the FSD data and 



/RMs(t) = fto + Ai sm I — 1 -I- Bi cos I — 
+A2sm( — j +B2C0S( — 



+e "cbo Acbo sin 



27rt 

-^cbo 



Brhn COS 



27rt 

^cbo 



, (32) 



where Oq is the average PSD profile width. The time bin- 
width was set to the cyclotron period, 149.185 ns, in order 
to average out oscillations in the time distribution caused 
by the bunched structure of the beam in the storage ring. 

The parameters Tcbo, Tcbo and T, the muon spin preces- 
sion period (21: /uj) are all fixed. In particular, T is set to 
the nominal value 4365.4 ns. The rms width oscillations 
contain precession frequency terms Ai and double 
precession frequency terms A2 and B2 due to the changes 
in average energy, detector acceptance and time-of-fiight 
of the decay electron within the precession period. CBO 
frequency terms A^ho and i?cbo arise because of variation 
in the decay electron time of flight with radius. 

Next, the vertical mean position versus time is flt to a 
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TABLE IV: Number of decay electrons used for PSD data analysis (in millions). 



Data Set 


PSD 1 


PSD 2 


PSD 3 


PSD 4 


PSD 5 


PSD (5 stations) 


Before B Correction n 


= 0.122 


66 


56 


61 


42 


70 


295 


Before B Correction n 


= 0.142 


61 


54 


60 


44 


71 


290 


After B Correction n 


= 0.122 


83 


73 


76 


65 


93 


390 


Total 


185 


157 


173 


139 


209 


975 



20^ 

E : 

E 

5 15- 
E, - 
a. 10- 



Chi2/ndf = 5.1 /4 
pO =1.04 ±3.34 




-5- 



(a)Before Br Correction 



25: 



£. 20: 



Chi2/nclf = 


3.6/4 


pO = 4.2 


± 5.4 




(b)Aftcr Br Correction 



FIG. 14: PSD analysis: EDM-sensitive parameter pi (see equation I31|l versus PSD station. Horizontal line indicates fit to a 
constant, pi is consistent with zero for both the data taken before and after Br was changed (see text for discussion). 



similar function: 



+A2g2 sin 



fAve{t) =yo+ Ag2 sin {^^Y^ + ^92 COS 

/47r<\ 

— ) + B2g2 COS 



2-Kt 

~Y 

Ant 



/ 27ri 

+e "cbo X 5cbo sm — h $cbo 

V-'cbo 

+Ccbo cos — h $cbo 

V Jcbo 



(33) 



where yo describes any misalignment between the beam 
and the vertical center of the PSD. There are precession 
frequency sine and cosine terms, Ag2 and Bg2 , and double 
precession frequency terms ^2^2 and i?2g2 all with the g — 
2 period fixed. The absolute phase of each detector was 
chosen so that the g — 2 number oscillation is described 
by Bg2 and any EDM signal would appear in Ag2 ■ There 
are also CBO sine and cosine terms, Acho and -Bcbo- The 
CBO decay time, period and phase are fixed by the RMS 
width fit. 

As in the FSD analysis, the very large false EDM signal 
caused by detector misalignment [2^ must be corrected. 
The EDM-sensitive parameter, Ag2, was plotted versus 
the PSD measured vertical detector offset 7/0 • The value 
of Ag2 at 2/0 = 0, obtained from a fit to a straight line, 
corresponds, once again, to the EDM-related amplitude 
that would be measured by a perfectly aligned detector. 



Fig. [T^ shows the Ag2 amplitude for the 5 PSD detectors, 
with separate points plotted for the data sets taken before 
and after the radial field correction. 



§. 20 

CM 

< 15 
10 

5 

-5 
-10 
-15 
-20 



Chi2 / ndf = 6.1/8 
pO =-0.75 ±1.68 
pi =-3.65 ±0.88 




Yo [mm] 

FIG. 15: PSD analysis: Fitted sine wave amplitude, Ag2 for 
the mean vertical position fit versus measured detector offset. 
A straight-line fit is overlaid. Since the intercept is consistent 
with 0, Ag2 is directly proportional to the measured detector 
offset. 
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B. Systematic Errors 

Systematic errors must be assessed separately for the 
two PSD analyses but most are similar to those of the 
FSD analysis. One of the largest systematic errors for the 
phase method is associated with the vertical alignment 
of the PSD with respect to the stored beam. To estimate 
the error, the relationship between the size of the false 
EDM signal and 1 mm of detector misalignment was es- 
tablished, using Monte Carlo simulation. The accuracy 
of the mean vertical position measurement was studied 
by imposing on the data a 10 % tile inefficiency and 5% 
gain variation for top/bottom halves of the calorimeter. 
The estimated error in vertical position is 0.2 mm. The 
Monte Carlo simulation was then used to determine the 
size of the false EDM signal produced per unit misalign- 
ment (see [1^), 140 /irad/mm. Multiplying that scaling 
factor by the estimated vertical misalignment yields an 
error of 29 /xrad on pi . The resulting error on the EDM 
is 0.74 X lO"^'' e-cm. 

Since the intercept of the sine wave amplitude ver- 
sus detector offset plot was used for the vertical position 
EDM measurement, detector misalignment was less im- 
portant to the overall systematic error. In that case, the 
systematic uncertainty was estimated using the slope of 
the sine wave Ag2 versus detector offset Yo plot combined 
with the error in the detector offset measurement. The 
systematic error for the detector misalignment is 1.7 /im 
in Ag2, resulting in an EDM error of 0.2 x 10^^^ e-cm, 
as in Eqs. [26] and [27l 

If the PSD is tilted with respect to the horizontal plane 
of the storage ring, horizontal beam oscillations produce 
apparent vertical oscillations in the average position mea- 
sured by the PSD. The PSD tilt was estimated from a 
two-dimensional plot of PSD Y v. X coordinates, for 
data taken more than 32 /isec after injection. The mean 
vertical position was calculated for each X tile and the 
ensemble of means was fitted vs. X to a linear function. 
The tilt of five PSDs was consistent with and no larger 
than 0.75 mrad on average. 

PSD tilt could also lead to asymmetric vertical losses 
which would affect the phase analysis. Higher energy 
electrons strike the calorimeter at larger radius, that is, 
closer to the storage region. Asymmetric losses could 
change the energy spectra and hence the phase of the 
detected electrons vs. vertical coordinate. The size of the 
error was estimated with the full Monte Carlo simulation. 
The induced EDM signal in parameter pi per mrad of 
detector tilt is 26 /xrad/mm (see '25^). Multiplying that 
scaling factor by the maximum tilt (0.75 mrad) yields an 
error 20 /irad/mm on pi. The associated error on the 
EDM is 0.5 X 10-19 e-cm. 

As noted before, a time-dependent energy scale differ- 
ence between the top and bottom of the calorimeter can 
induce a false EDM signal. The energy balance procedure 
was checked by examining events in which the signal was 
restricted to the top or bottom of the calorimeter. The 
measured average top/bottom energy scale difference is 



less than 0.1% after the time-dependent energy correction 
was applied. Monte Carlo simulation indicates that a 1 
% top/bottom gain imbalance produces a 43 /xrad/mm 
error, corresponding to an error of 4.3 /irad/mm in our 
case. 

In order to study the error due to tile inefficiency, an 
artificial 10% inefficiency for tile 3 (for all PSDs) was 
imposed on the data. (Tile 3 corresponds to the upper- 
most tile used in the analysis). The intercept change in 
Fig. [iniis only 1.8 nm, which is taken as an estimate of 
the systematic error for the vertical mean method. 

Errors for the phase method are summarized in Ta- 
ble |Vl The total uncertainty is 50 /irad/mm on pi, 
1.3 X 10^19 e-cm on the EDM. Systematic uncertain- 
ties for the mean vertical position vs. time method 
are presented in Table IVII The total systematic un- 
certainty for the vertical position method is 9.6 /im (or 
1.1 X lO-"'^^ e-cm) and statistical error is 1.7 //m (or 
0.2 X 10^19 e-cm). 

C. PSD Summary 

Results from the two PSD analysis methods are con- 
sistent with zero and with each other. Averaged over the 
five PSD stations, the EDM signal from the phase anal- 
ysis is = (-0.48 ± 1.3) X lO^^^ e-cm, where 0.73 x 
lO^i^e-cm is the statistical error and 1.09 x 10"^^ e-cm 
is the systematic error. Averaged over the five PSDs, 
fitting the PSD vertical position versus time yields an 
EDM value of d^- = (-0.1 ± 0.73) x lO'^^ e-cm, where 
0.28 X 10-19 e-cm and 0.70 x lO-^^ e-cm are the statistical 
and systematic errors, respectively. 

VII. COMBINATION OF EDM RESULTS 

The traceback data set is entirely independent of those 
used for the FSD and PSD analyses and there are min- 
imal correlations between the systematic errors of the 
traceback analysis and those of the other two. 

The combined measurement from the traceback and 
FSD analyses on the /i+ is = (-0.1 ± 1.0) x lO-^^ e- 
cm, providing a limit on the permanent electric dipole 
moment of the positive muon of 

\df,+ I < 2.1 X 10-19 e - cm (95% C.L.). (34) 

The /*- result is taken from the PSD mean vertical po- 
sition analysis: d^- = (—0.1 ± 0.7) x lO-^^ g.cm with 
corresponding limit 

I < 1.5 X 10^19 g . cm (95% C.L.). (35) 

It is important to note that the errors on the traceback 
measurement are entirely statistical while those of the 
FSD and PSD measurements are dominantly systematic. 
Although the increased vertical segmentation of the BNL 
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TABLE V: Errors for the Phase Method. The first six entries of the central column list the size of a physical error and the 
sensitivity (assumed to be linear) of the EDM result to that error. The column to the right gives errors on EDM-sensitive 
parameter pi and in the totals, the error on the EDM itself. 



Source 


Sensitivity 


Result 


Detector Tilt 
Detector Misalignment 
Energy Calibration 
Muon Vertical Spin 
Radial B field 
Timing 

Total systematic 
Total statistical 


26 ^rad/mm/mradx 0.75 mrad 
138 ^rad/mm/ mm x 0.2 mm 
43 ^irad/mm/ % x 0.1% 
1.0 /irad/mm x 8% 
0.72 /Ltrad/mm/ppm x 20.0 ppm 
17.0 /Ltrad/mm/ns x 0.2 ns 


20 fi rad/mm 
28 /i rad/mm 

4.3 /i rad/mm 
8.0 jj. rad/mm 
14.4 jj, rad/mm 

3.4 /i rad/mm 

38 ^irad/mm (0.93 x 10"^'' e-cm ) 
28 ytirad/mm (0.73 x 10"^^ e-cm ) 


Total 




47 //rad/mm (1.2 x 10"^^ e-cm ) 



TABLE VI: Errors for the Mean Position Method. The first five entries of the central column list the size of a physical error 
and the sensitivity (assumed to be linear) on the residual value of the EDM-sensitive parameter Ag2 for a nomiminal detector 
offset of 0. In the last three rows, the systematic, statistical and total errors on the EDM axe also given. 



Source 


Sensitivity 


Result 


Detector Misal. 


4.1 /im/mm x 0.4 mm 


1.6 fim 


Detector Tilt 


3.0 /xm/mrad x 0.75 mrad 


2.3 /xm 


Energy Calibration 


28.0 urn/ %x 0.1 % 


2.8 urn 


Muon Vertical Spin 


0.4 fim/% X 8% 


3.2 fim 


Radial B field 


0.13 /jm/ppm X 20 ppm 


2.6 fim 


Timing 




1.5 fim 


Tile ineff. 




1.7 /ttm 


Fitting 




1.0 jim 


Total Systematic 




6.2 nm (0.7 X 10"^^ e-cm ) 


Total Statistical 




1.7 urn (0.2 X IQ-i^ e-cm) 


Total 




6.4 /im (0.73 x 10"^^ e-cm) 



detectors is a distinct improvement over the two-paddle 
CERN III approach, significant further progress with this 
basic technique would be exceedingly difScult. By con- 
trast, relatively simple remedies - increasing the geomet- 
ric acceptance and reducing or eliminating the recovery 
period after injection - would lead to direct improvements 
in the traceback result. 

In determining joint results from the /i+ and fi^ data 
sets, we must consider the correlation of errors in the 
FSD(/i+) and PSD(/x-) analyses. While the FSD and 
PSD data sets arc entirely independent, three of the sys- 
tematic errors: the vertical spin, top/bottom calibration 
and radial B-field errors are correlated and, although it 
is impossible to judge the extent of the correlation, one 
might assume that the B- and E-fields of all bending 
and focusing elements are exactly reversed when chang- 
ing from positive to negative beam and that these three 
errors are fully correlated. 

The EDMs of the /i+ and /it" are found to be in ac- 
cord with the CPT requirement that = — rf^j- , or, in 
terms of directly measured quantities, 77^+ = rj^- . The 



result is 

d^,- + = (-0.2 ± 1.3) x 10"^^ e • cm. (36) 

In summing the results for //+ and /x~ to construct this 
difference, the correlated systematic errors in the FSD 
and PSD analyses should tend to cancel, as described 
above. However, because of uncertainties in the extent 
of the correlation, as a conservative measure, we take 
them to be completely uncorrelated. 

Under the assumption of CPT invariance, we take a 
weighted average of the result and the opposite of 
the fi'^ results. In this case, once again as a conservative 
measure, we take the three correlated errors of the FSD 
and PSD analyses to contribute maximally to the total 
error. All other errors are added in quadrature. We 
obtain 

rf^ = (-0.1 ±0.9) X 10"^^ e-cm. (37) 
This corresponds to the limit 

\d^\ < 1.9 X 10^1^ e - cm (95% C.L.), (38) 
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approximately a factor of 5 improvement over the previ- 
ous limit. A summary of the results from each detector 
system, as well those from CERN III, is shown in Table 

Em 

The observed difference between the standard model 
and the experimentally determined values of the muon 
anomaly, Ja^ = 295(88) x 10"" could be at- 

tributed, in principle, to a shift in the muon spin preces- 
sion frequency caused by a non-zero EDM. See Eq. p6)) . 
This would lead to an EDM of dsa^ = ±2.39(0.36) x 
10"^^ e-cm. The new limit on the EDM (Eq. [37l) would 
predict (5a^(EDM) = 0.5(-H0.51/ - 0.33) x 10"". The 
probability that d^Q^ is compatible with the limit from 
Eq. [37] is 2%, suggesting it is unlikely that fo^ arises 
from a non-zero EDM. 
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TABLE VII: Summary of EDM results from the CERN III experiment and from the three analyses of E821. The units are 
xlO~^®e- cm in all cases. 
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